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Abstract 

The photoproduction reaction 7^ — > fi + n~p has been studied in ep interactions 
using the ZEUS detector at HERA. The data sample corresponds to an integrated 
luminosity of 43.2 pb _1 . The T meson has been observed in photoproduction 
for the first time. The sum of the products of the elastic T(15), T(2S), T(35) 
photoproduction cross sections with their respective branching ratios is determined 
to be 13.3 ± 6.0(stat.)~^23( s y s t-) pb at a mean photon-proton centre of mass energy 
of 120 GeV. The cross section is above the prediction of a perturbative QCD model. 
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1 Introduction 



Perturbative QCD can be applied to ep scattering to calculate the amplitude for elastic 
production of heavy vector mesons. Previous HERA results on J/ip meson production 
at Q 2 ~ and for 0.25 < Q 2 < 40 GeV 2 |], H nave shown that the rise of the cross 
section with W as well as the dependence of the cross section with Q 2 can be described 
by perturbative QCD models || £|, |J. Only upper limits of a ■ B exist in the literature 
for inclusive T production in lepton-hadron collisions 0. 

In this paper we extend the study of elastic photoproduction of vector mesons at HERA 
to T mesons. The measurement is made with the ZEUS detector using a data sample 
which corresponds to an integrated ep luminosity of 43.2 pb _1 . The improved luminosity 
allows the study of the reaction 7p — ► fi + fi^p for invariant masses beyond the T 

mass region. A first measurement of cr 7P ^T P • B in photoproduction is presented in the 
kinematic range of the photon-proton centre of mass energy 80 < W < 160 GeV, where 
B is the T branching ratio to muons. The T(15), T(2S) and T(3S) resonances are not 
resolved in this measurement. Under the assumption that the production ratios of Y(IS'), 
T(25) and T(35) are the same as those measured in hadron-hadron collisions J5], ||, we 
determine the T(IS') and the ratio of the T(15) to J/ip photoproduction cross sections 
and compare them to a pQCD inspired model ||. 



2 Experimental Conditions 

In the years 1995-97 HERA collided positrons of 27.5 GeV with protons of 820 GeV, cor- 
responding to a centre of mass energy y/s = 300 GeV. A description of the ZEUS detector 
can be found in references || W\. The primary components used in this analysis were the 
central tracking detector, the uranium-scintillator calorimeter, the muon chambers and 
the proton remnant tagger. The central tracking detector (CTD) 1TIJ operates in a 1.43 T 
solenoidal magnetic field. It is a drift chamber consisting of 72 cylindrical layers, organized 
into 9 superlayers. It was used to identify the vertex and to measure the momenta and di- 
rections of the muons. The transverse momentum resolution is u{pt) jpt — [0.005p t ] ©0.016, 
with p t in GeV, for full length tracks. The calorimeter (CAL) ]12| is hermetic. Under 
test beam conditions, it has energy resolutions of 18%/\/E for electrons and 35%/\/E for 
hadrons. The time resolution is below 1 ns for energy deposits greater than 4.5 GeV. The 
CAL was used to reject cosmic rays and beam halo background by timing and to identify 
minimum ionizing particles (m.i.p.). It is surrounded by a magnetized iron yoke with a 
field of 1.4 T produced by conventional warm coils. The muon system (MUO) consists of 
tracking detectors (forward, barrel and rear muon chambers: FMU HlQ"l , B/RMU [ |i~3|j ) 
placed inside and outside the yoke covering the polar angles from 10° to 171°. They were 
used for the trigger and to tag the muons by matching segments in the MUO chambers 
with tracks in the CTD and m.i.p. 's in the CAL. The proton remnant tagger (PRT) con- 
sists of two stations of scintillation counters surrounding the beamline at Z = 5 m and 
Z = 24 m [], and tags particles in the forward proton direction in an angular range 6 to 26 

lr rhe right-handed ZEUS coordinate system is centred on the nominal interaction point (Z — 0) and 
defined with the Z axis pointing in the proton beam direction, and the horizontal X axis pointing towards 
the centre of HERA. 
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mrad and 1.5 to 8 mrad, respectively. It was used to estimate proton dissociative back- 
ground. The luminosity was measured to a precision of 1.5% from the rate of energetic 
bremsstrahlung photons produced in the process e + p — > e + jp. 

3 Kinematics 

Figure |l|a shows a schematic diagram for the reaction 

e+(k) p{p) -> e + (k') V(v) p(p'), (1) 

where the symbol in parenthesis denotes the four-momentum of the corresponding particle, 
and V indicates a J /if) or an T. 

The kinematics of the inclusive scattering of unpolarised positrons and protons are de- 
scribed by the positron-proton centre of mass energy squared (s) and any two of the 
following variables: 

• Q 2 = —q 2 = — (k — k') 2 , the negative four-momentum squared of the exchanged 
photon; 

• U = (Q ' P)/(k ' P)i the fraction of the positron energy transferred to the hadronic 
final state in the rest frame of the initial state proton; 

• W 2 — (q + p) 2 = —Q 2 + 2y(k ■ p) + M 2 , the centre of mass energy squared of the 
photon-proton system, where M p is the proton mass. 

For a complete description of the exclusive reaction e + p — > e + V p (V — > ft + ft~) the 
following additional variables are required: 

• My, the invariant mass of the pair; 

• t = {p — P') 2 ) the four-momentum transfer squared at the proton vertex; 

• $, the angle between the vector meson production plane and the positron scattering 
plane in the photon-proton centre of mass frame; 

• 6 h and (fi h , the polar and azimuthal angles of the positively charged decay lepton in 
the V helicity frame. 

In this analysis, photoproduction events are selected by requesting a pair from the 

interaction point and nothing else in either the CTD or the CAL. For the selected events 
the Q 2 value ranges from the kinematic minimum Q^ in = M 2 y 2 /(1 — y) m 10~ 9 GeV 2 , 
where M e is the electron mass, to the value at which the scattered positron starts to be 
observed in the uranium calorimeter Q 2 max ~ 1 GeV 2 , with a median Q 2 of approximately 
5 x 10~ 5 GeV 2 . Since the typical Q 2 is small, the photon-proton centre of mass energy 
can be expressed as 

W 2 w AE p E e y w 2E P {E - p z ) v , (2) 

where E p is the laboratory energy of the incoming proton and (E —pz)v is the difference 
between the energy and the longitudinal momentum of the vector meson, V, as determined 
from the CTD tracks assuming the muon mass. $ is not measurable since the scattered 
positron is not detected. 
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4 Event Selection 



Elastic events were selected using dedicated triggers. Trigger cuts are superseded 

by the following offline cuts: 



• CAL timing and reconstructed interaction vertex consistent with the nominal ep 
interaction to reject non e + p background; 

• two oppositely charged tracks from the vertex, at least one of which matches a 
segment in the B/RMU chambers or a hit in the FMU chambers, and no other 
track in the CTD; 

• at least 3 CTD superlayers per track, limiting the polar angular region from ~ 17° 
to ~ 163°; 

• acollinearity of the tracks Q < 174° where Q is the angle between the two tracks, to 
reject cosmic ray events; 

• invariant mass of the two CTD tracks, treated as muons, larger than 2 GeV; 

• CAL energy associated to each track consistent with the energy deposit of a m.i.p.; 

• no CAL cell (apart from those associated with a /i) with energy greater than 150 
MeV, well above the CAL uranium noise level. 



The events were then selected in the kinematical range of W between 80 and 160 GeV, 
corresponding to an acceptance of ~ 40% for the T; the acceptance for the J/ip falls 
smoothly from 38% to 10% with increasing W over this W range. 



5 Monte Carlo Simulation 



To compute the acceptance, the reaction [I] (Figure |l|a) was simulated using the DIPSI 
generator [14|. DIPSI is based on a model developed by Ryskin 0] in which the exchanged 



photon fluctuates into a qq pair which interacts with a gluon ladder emitted by the incident 
proton. The parameters of the model are the strong coupling constant a s (assumed fixed) , 
the two-gluon form factor and the gluon momentum density of the proton. The cross 
section dependence on W and t is fixed by these parameters. 

Proton dissociative vector meson production, e + p — > e + VN (Figure |I|b), was simulated 
using the EPSOFT [I5| and the DIFFVM [TO generators. Both are based on the assump- 
tion that the diffractive cross section at large is of the form d 2 a / d\t\dM% oc e~ b '*'/M^ 
where is the mass of the dissociative system; we have used 6 = 1 GeV~ 2 and (3 = 2.2. 
The simulation of the dissociative system includes a parametrisation of the spectrum in 
the resonance region which differs for the two generators. 

The background Bethe-Heitler process, in which a lepton pair is produced by the fusion 
of a photon radiated by the positron with a photon radiated by the proton, was simu- 
lated using the LPAIR generator |T7[ . Both elastic and proton dissociative events were 
generated. 
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All Monte Carlo events were passed through a simulation of the ZEUS detector and trigger 
based on the GEANT program |18j and analysed with the same reconstruction and offline 
selection procedures as the data. The overall acceptance in a selected kinematic range 
was obtained as the ratio of the number of Monte Carlo events passing the cuts to the 
number of events generated in the same range. The acceptance, calculated in this man- 
ner, accounts for the geometric acceptance, for the detector, trigger and reconstruction 
efficiencies, and for the detector resolution. 



6 Analysis 

The measured mass distribution, for the sample of events obtained by the selection 

described in Section [|, is shown in Figure 0. The signals of the J/tp, i>(2S) and T's 
(unresolved) are apparent. The continuum is well described by the Bethe-Heitler process 
apart from an enhancement in the region of 6 GeV, which is consistent with being a 
fluctuation. 



6.1 Extraction of the Signal 

The elastic Bethe-Heitler process can be calculated to good precision since it is a pure 
QED reaction. There is, however, some arbitrariness in the parametrisation of the proton 
dissociative component. Comparing the data with the LPAIR simulation in the mass 
interval 4 to 8 GeV, where we expect the Bethe-Heitler process to dominate, we find 
that the simulation underestimates the energy deposited in the forward region of CAL 
due to proton dissociation. This leads to an overestimate of the Bethe-Heitler events 
which pass the selection cuts if both the elastic and the proton dissociative processes 



are normalised to luminosity by their calculated cross sections [|17J . However, since the 
shape of the mass distribution predicted by LPAIR is the same for both processes, to 
evaluate this background we have normalised Monte Carlo and data in a mass window not 
containing resonances (4.2 to 8.4 GeV). This corresponds to adding a 10% contribution to 
the elastic Bethe-Heitler distribution normalised to the luminosity. In Figure |2] the data 
are compared with this renormalised Bethe-Heitler distribution. The spectrum outside 
the resonance regions is well reproduced and the distribution has been used to subtract 
the background under the resonances. 

The limited statistics and the ^ 4 ' \x~ mass resolution of 0.3 GeV in the T region do not 
allow to distinguish between the T(15), T(25) and T(35) states. The mass window 8.9 
to 10.9 GeV (i.e. from twice the resolution below the T(1S) nominal mass, 9.46 GeV, 
to twice the resolution above the T(3S) mass, 10.36 GeV fl9fl ) was chosen to count the 
T events. In this region there are 57 events while 39.9 are estimated to be background. 
The mass spectrum in the T region after background subtraction is shown in the insert 
of Figure the mean mass and the r.m.s. for the background corrected signal in the 
window are 9.9 ± 0.2 GeV and 0.47 GeV, respectively. The simulation of the detector 
response to a mixture of T states J7| yields a mean mass of 9.7 GeV and an r.m.s. of 0.42 
GeV. The r.m.s. is smaller since background was not considered. A total of 4257 events 
are counted in the J/ip mass window from 2.8 to 3.35 GeV while 306 are expected from 
the Bethe-Heitler process. 
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The amount of proton dissociative resonant background remaining in the sample after 
the cuts described in Section f| is estimated using the number of events tagged by the 
PRT and the EPSOFT Monte Carlo to correct for the detector acceptance. In the J/if> 
mass region the fraction of dissociative events is determined to be 25 ± 2(stat.)% . For 
T production, the same fraction of proton dissociative events has been assumed. Within 
the limited statistics, this assumption is consistent with the number of events in the T 
mass region tagged by the PRT. 

6.2 Systematic Uncertainties 

A study of the systematic uncertainties on the measurements has been performed. The 
systematic uncertainties, listed below, have been divided into two classes, those common 
to the T and J/ip analysis, which cancel in the measurement of the ratio, and those 
remaining. 

• Common systematic uncertainties: 

— uncertainty on the CTD first level trigger efficiency, ±5%; 

— uncertainty on the muon chamber and muon trigger efficiency, ±10%; 

— uncertainty on the proton dissociation background estimated using DIFFVM, 
(+9,-2)%. 

• Systematic uncertainties specific to each meson: 

— using a different mass region for normalisation of background contribution, 
from 4 to 6 GeV, +10%; 

— varying the T mass window by ±300 MeV, ±10%; 

— varying the cos 9h distribution used in the Monte Carlo from that required by 
s-channel helicity conservation (1 + cos 2 9^) to a flat distribution for the T and 
between the limits allowed by the existing measurements for the J/ip, —8% for 
the T, -4% for the J/ip; 

— the uncertainties on the J/ip and T muonic branching ratios, ±3%. 

All of these uncertainties are added in quadrature yielding an overall systematic error on 
cr • B of (+20, —17)%. The systematics which affect the a 1P ^ p / a 1P ^j/^ p ratio are only 
those included in the second group which give contributions of (+14, —13)%. 

Other possible sources of systematic uncertainty (variation of the energy thresholds used 
to select the elastic events, removal of the m.i.p. requirement on the energy associated with 
the tracks, variation of the J / ip mass window, uncertainty in the luminosity determination, 
variation of the DIPSI parameters, use of an unbinned fit to extract the T signal) give 
negligible contributions. 
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7 Results and comparison with QCD predictions 



We have calculated a ■ B as 

_ N evt (1 - T ) 
a e +p^ e + V p ■ o — ~~X~£ ' ^ ' 

where N evt denotes the background subtracted number of V signal events, T is the esti- 
mated fraction of proton dissociative events, A the acceptance, £ the integrated luminosity 
and B the muonic branching fraction. Photoproduction cross sections were determined 
by dividing the electroproduction cross sections by the photon flux, calculated according 
to reference |20| . 

The production cross sections multiplied by the muonic branching ratios for the unresolved 
T(IS'), T(25) and T(3S') are summarized in Table fl]. cr 7p ^ Tp • B is determined to be 
13.3 ± G.O^stat.^l'l^syst.) pb at a mean photon-proton centre of mass energy (W) = 120 
GeV. The first error is statistical and the second is the systematic uncertainty. In Table [I] 
we report also the ratio a 1P ^rp • B summed on the T states over cr 7P ^j/^ p in the W range 
between 80 and 160 GeV. 

The QCD-based model of reference H gives predictions for the hard diffractive photo- and 
electroproduction of heavy vector meson J/ip and T(IS') within the leading logarithmic 
approximation. To compare the data with the predictions, the experimental cross section 
for T(15') has to be derived. As our data does not allow the relative fractions of T(1S), 
T(2S) and T(3S) to be determined, we have to make assumptions on their relative pro- 
duction rates. Assuming that the cross sections times branching ratios are the same as 
measured by CDF §, a ■ B(T(2S))/a ■ B(T(1S)) = 0.281 ± 0M0(stat.) ± 0M8(sys) and 
a-B{T(3S))/a-B{T{lS)) = 0.155 ±0.024(stat) ±0.021(sys), the T(15) accounts for 70% 
of the signal. Using the muonic branching ratio B(T(1S) — > = (2.48 ± 0.007)% 

|fT9[ 1, we derive 



a w ^ r{1S) p = 375 ± 170(stat.)± 7 Jl(syst.) pb at W = 120 GeV 

and the ratio 

^ 7 p^t(is) p/v 7 p^JW p = (4.8 ± 2.2(stat.)tol(syst.)) x 10" 3 . 

In Figure ^|(a),(b) these values are compared with the theoretical calculations which give 
~ 60 pb for <t 7P ^t(is) P and ~ 0.001 for the ratio, both weakly dependent on the structure 
function parametrisation used. Our measurement is higher than the predictions. 
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W range 
(GeV) 


(W) 
(GeV) 


N evt 


A 


(Pb) 




(Pb) 




80 - 160 


120 


17.1 ± 7.5 


0.43 


0.68 ± 0.30±g;i| 


0.051 


13.3 ± 6.0±|;| 


1.7 ± 0.8±g;| 



Table 1: The results for the unresolved T cross sections times the muonic branching ratios, 
Y^i Ci ■ where % runs over the three states T(IS'), T(2S), T(3S) and £>j is the branching ratio 
of each state to muons. N evt is the number of events after subtraction of the non-resonant 
background contribution: 75% of N evt is attributed to the elastic reaction, the remaining 25% 
to the proton dissociative process. A is the acceptance and $t is the photon flux used to 
calculate the 7J9 cross section (J 1P ^ P from the ep cross section a ep ^ e r P - The last column 
contains the ratio of unresolved T cross section times the muonic branching fraction to J/ip 
cross section. The first uncertainties are statistical and the second ones systematic. 
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Figure 1: Schematic diagrams for diffractive vector meson (V) electroproduction. (a) 
Elastic V production, (b) Proton dissociative V production where the proton dissociates 
into a hadronic system of invariant mass M^. The dotted line represents a colorless 
exchange between the virtual photon and the proton. 
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Figure 2: Mass distribution of pairs. The histogram represents the simulated 

Bethe-Heitler background. Points in the J/ip region are connected by a dotted line to 
guide the eye. The insert shows the signal remaining in the T region after subtraction of 
the non-resonant background. 
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Figure 3: (a) 0" 7P _^x(i5) P and (b) the ratio cr 7P _>Y(is) p/^-yp-*j/<s p as a function of W, 
the centre of mass energy of the photon-proton system. Data (full squares) are compared 
with the predictions of for the GRV94(HO) (dashed line) and CTEQ(4M) (full line) 
parametrisations of the proton structure function. The data have been scaled taking into 
account the muonic branching ratio and the contributions of the T(2S) and T(3S), as 
described in the text. The inner error bars show the statistical errors, the outer bars 
correspond to the statistical and systematic errors added in quadrature. 
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